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ABSTRACT 

A computer-assisted course was designed to provide 
students with an understanding of modeling and simulation techniques 
in guan titiative ecology. It deals with continuous systems and has 
two segments*. One develops mathematical and computer tools, beginning 
with abstract systems and their relation to physical systems. 
Modeling principles are next presented, with emphasis on 
compartment-flow models, followed by discussion of simulation. The 
FOBTRAN IV language and the IBM 360 Continous System Model Program 
are then described in detail, and equilibria and stability ara 
treated. In part two, students were burdened by the math required and 
the amount of time spent on programing languages. The first 
difficulty was overcome by the introduction of cpmpartment-flow 
models, the second by the development' of an easy-to-use computer 
program. This Compartment Flow Model Simulation Program (COFLO) is an 
input language which takes a description of a cqmpa rtment-f low model 
and simulates the behavior of the modelled system, converting the 
model system into a FORTRAN program for compilation and execution. 
COFLO has been implemented on the UNIVAC 1106 computer at the 
Southeast Regional Data Center of Florida Atlantic and Florida 
International Universities. (LB) 
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COFLO: 



A ConPUTHR AID FOR TEACHIKG ECOlor,lCAL LiIJiULAlJON 



Dopilta'^nt of Katl.<?iriatirs 
Florida Atlantic University 
Hoca luton, Florida 32^^2 
(J05) 395-5100 



Ir this papor the author de£crib*=r3 hl:5 exp^-rivnces t<^achinq a course- in ffiod»^Iin<; and 
simulition methods in ^colcqy Florida Atlantic University durir.9 th€ rprir.g quart^^r ol tU*- 
1971-1^72 academic year. H-? also describes a n<^w i:iDulat icr; lan^juatj'? c*?si<jn^d F.i)»jci t ica 1 1 y 
for stident uKe in ecological simulation which was d€Vt>lop#;d ti.i an cutcjrowth of tht cours^?. 



The course was scaievha*. unu£ual in that it brought tog^th'^r a <jroup of i^tuden*.s with 
diverse backgrounds, — natheniat xcs, computer science, cr hiolcgy-*and attempted to give rheta lu 
understanding of those modeling and si mala t ion techniqu-?:: which are r\Oi:t rel-vart to 
quantitative ecology. The need for such a courre i:i obviour on iiov-^rdl grounds. (a) 
Gystens o ""ogy is an ea.ergent discipline whoso full dev-i lo{:3ienr will retjuir*^ con trilnit ions 
from (oa theiijw cicia ns and computer ccieiitijts as well as ecolcgi:>t3. (b) Rapid developnent of 
systeas ^colo.jy must be encouraged to facilitate Jc;.<; range planning to protect the 
environoent. (c) The problems of systeas ecology ar^ ^uft iciertly different troni those 
encountered in modeling and simulation in busine^.'^ and engineering tliat a separate couri;e is 
required for modeling and simulation methods in ecology. 

The course was concerned pricarily with codcling ccntinuouc Gynt<?:^r; vith only flight 
attention to discrete and probabilistic systems. It was divided into 'wo s'oc.jnts. The 
first dealt primarily with the development of the mathematical and computer tools required 
for continuous system simulation, and the second segment ccnsist^rd of an examination of 
several papers in systems ecology to permit observation cl applications 01 the tools which 
had been developc-d. The principal source of material for t te first segment of th<:- course was 
part I of Systems Anal^ris and Simulation in Ecology, volume I, edited by Bernard C. Patten 
[3]. Papers for the second segment of the course w€re selected primarily fio:n the remainder 
of Patten and from Heinmets f 1 ). 

The first segment essentially followed the form of part I of Patten [3]. The notion of 
an abstract system was introduced, and the features of such systems and th^^ir relation to 
physical systems was discussed. The principles of modeling were discussed next with eraphasis 
on cOBpar t Dent ^ flow models. This having been done, the notion of simulation was introduced. 
First, analog computation was discussed, and then digital cocf utatiofi. Since the class did 
not have access to an analog computer, the IBM 1130 continuous System Modeling Program, which 
provides a digital computer simulation of an analog computer, based on input containing a 
description of the analog computer circuit, was ustd. The FOFTRAN IV programming language 
was briefly described, and then the IBM SysteBi/360 Continuous System Modeling Program was 
presented in detail. The first segment of the course concludod with a brief dii;cuGsion of 
equilibria and stability. 

In the second segment each student was given the task of pjoviding a bri<^f oral 
presentation of a published paper in systems ecology for the rest of the clas:^. 

The most striking impression of the material of the course was that the non- 
mathematicians were heavily burdened by the ma themat ics- -part icula r ly differential 
equations — and that an extremely large amount of tim-? van spent on urogrartming languages at 
the expense of other material. Further analysis of the mdtetidl in the coucs-; revealed ttiat 
the use of compartment-flow models, which are quite natural for many ecological systerr*3, 
could permit a reasonably full development of sir?iilatior concept;, without th^ explicit 
introduction of differential equations, provided that a mechanism was available for the 
direct simulation of a model expressed in ccmpar tmen t-£ low form. Fur th'^r mor the 
developffient of an easy to use computer program to handle this problem would substantially 
reduce the amount of time spent cn programming languages. 
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The Compa rtoent-Flow flod^l Sioulation Progran (COFLO) was designed to fleet the 
requircaent set forth above, while .illoving for the greatest possible variety of 
applications. COFLO is an input language which takes as its input a description of a 
coBpart Bent-flow model and sioulates the behavior of the oodeled systecj for a given tine 
period^ printing the results at prescribed time intervals. It functions by converting the 
Dodel description into a FORTRAN program which it then causes tc be compiled and executed. 

The program requires that the model have the following form. The non-conservative 
coopartsen* is labeled 0, and each of the other compartaents is labeled with soxe in'^.eger 
between 1 and 99. The coapartotiit labeled i is associated with the variable Xi. The 
independent variable may be referred to either as T or TIHE. 

There are five control statements: ♦NAMES, ♦FLOWS, ♦INITIALIZE, ♦TIMER, and tRUS. The 
control statements appear in the input in the order indicated above. The first four are each 
followed by model description statements. 

♦KAflES is used to indicate that the following statements describe naaie.j with a form 
other than Xi which are to be used for the variables associated with certain compartments. 
The statements have the form 

Xi = name 

or 

T = name. 

The ♦NAMES section is optional and is only present if non-standard names are used. 

♦FLOWS precedes the actual description of th<? raodc-1 in terms of the individual flow 
equations. The equation for the flow from compartment i to compartment j is entered in the 
form 

F(i/J) = expression. 

The expression is written using the usual FORTRAN conventions. 

♦INITIALIZE indicates that the following statements provide initial values for 
compertaent variables or values for computational constants. These statements have the 
form 

variable = expression. 

Any variable which is set otherwise ifiitialized is initialized tc zero, and if all variables 
are to be initialized tc zero the section may be omitted. 

The *TIMER section is used to set START and STOP values for the timer, and DELTA which 
detorainei; the interval at which results are printed. The values are entered in the form 

name = value, 

and if value is zero the entry may be omitted. 

ihe ♦RUN card indicates the end of the model description and starts the simulation run. 

COFLO has been implemented on the UKIVAC 1106 computer at th-j Southeast Regional Data 
Center of Florida Atlantic and Florida International Universities. It may be used either 
from cards in the batch mode or from the interactive terminals in demaiid mode. In the later 
case the user receives immediate diagnostics if errors are made in entering the model and 
receives the results of the simulation without the delays inherent in batch utilization. 

The operation of COFLO is bost illustrated by example. Ir. Appendix I we have a COFLO 
run foe a partially forced non-linear version of odum's food chain model for Silver Springs, 
Florida [2], [3]. The corapartmonts are numbered as fellows: 







2 - 


" herbivoros 


3 - 


- car ni vc::e 


U - 


- top carnivores 


5 - 


- d'^rCOQposors 



The paraneters Tij are feeding rates, Ri ar^r respiratior. rates, ai are natural nort lity 
rates, and Li ate Icsses ^lownstr^an. 

rte prograo may also be used for nodels i/bijh do not fit the ccopa rtment-f low model form 
as illustrated by the exaniple in Appendix II whxch uses the LotKa-Volterra equations. 

Several impro ven^^nic are plajined tor COFLO in the near futur*. In particular a plotin<j 
capability is planned and a library of special functions will be adoed. Copies of tho COFLO 
code and u^jors guide will be oade available in the near future. 
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APPENDIX I 



>XOT .COFIO • SllVc;^ SPRINGS KODEL 



COMPAKTKENI - FLOW SIMULATION' PROGRAM 



•FLOt;S 

F(?»2) =K? 

F( 1 -T1?*X1*X? 

F( 1 id)=«l*Xl 
F(2«5)=M?*X? 
F(3i?:)=K.>X3 
p(4i5>=M4«X4 
F( Ii J)=LR1*X1 

F(3f'J)=R.?»X3 

F(4»2)=R4»X4 

F(5»e)=R5*X5 

•rM riALtZE 

X1=1421«26 

X2=?13.44 

X3=62.?6 

X4=S.87 

X5-24.3e 

Fl=,?- 

F2=496. 

T12=2C74./(X1»X2) 

123=382. /(X2^X3) 

(34=^21. /(X3*X4) 

KU34G5./X1 

«2=1205./X2 

M3=46./X3 

f<4 = 6./X4 

Rl=n374./Xl 

LI =2408#/Xl 

PL1=R1*L1 

R2=199N/X2 

R3=31 7./X3 

R4=13./X4 

R5=4593*/X5 

•TIMER 

SI OP =13. 

0ELTA=3.5 

TAOLE OF IN'CLUDEO FLOWS 
1 2 3 4 5 
2 XX 

1 X X X 

2 X XX 

3 X XX 

4 X X 

5 X 



END OF PREPARATION PHASE 



EAOT 





• 00 


X{l) = 


.3421*04 


.2134+03 


.6206+ 


02 




•2438+02 




.50 




•1577^04 


•7233+02 


.2338+ 


02 


• 72nt Ji 


.1059+02 


1 


• 22 


x(i) = 


.8520^03 


.4735+02 


.2242+ 


01 


•2J03+01 


•4446+31 


1 


• 50 




.4720^03 


.4072+0? 


.2173+03 


.1013+01 


.3257+01 


2 


.00 


x(i) = 


•2638^03 


•3783+02 


•2031- 


01 


•3477+00 


. 1 860+1^ 


2 


.50 


x<i) = 


.143^^03 


.3640+02 


•1853- 


02 


.1192+00 


.1673+01 


3 


.00 


x(i) = 


.8323^02 


.3565+02 


•1669- 


03 


•4383-01 


.1779+00 


3 


.50 


X( 1 ) = 


•4CS5^02 


.3525+02 


•1493- 


04 


•1399-01 


.1164+01 


4 


.00 


x(i) = 


.2039^02 


.3502+02 


.1330- 


05 


•4705-02 


.5923 ♦•03 


4 


-50 


x(i) = 


.14S7f02 


•3490+02 


.1182- 


26 


•1643-02 


.1010+01 


5 


.00 


x(i)=^ 


.8378^01 


•3433+02 


«1049- 


\dl 


•5629-03 


.9019+00 


5 


.50 


X(!) = 


.'.722 + 01 


•3479+02 


•9300- 


09 


•1929-03 


.9718+00 


b 


.00 


X(l) = 


•2661^01 


.3477+32 


•8244- 


10 


.5610-04 


.9533+00 


G 


♦ 50 


x(i) = 


.1503+01 


.3476+02 


•730C- 


1 1 


.2265-04 


.9535^00 


7 


♦ 00 


X( 1 ) = 


• 3455 + 00 


.3475+02 


•6474- 


12 


• 7761 


.9 1 91+32 


7 


.50 


x(i) = 


.4765+00 


•3474+02 


•5736- 


13 


• 2559-1^5 


.6187+00 


8 


.00 


x(i) = 


•2686+03 


•3474+02 


•5382- 


14 


.9112-06 


.9340+00 


8 


.50 


x(i) = 


•1514+00 


•3474+02 


•4503- 


15 


.3122-06 


.9454+00 


3 


.00 


x(i) = 


•3535-01 


.3474+02 


•3939- 


16 


•1070-06 


.0455+00 


9 


.50 


X(l) = 


•4811-01 


.3474+02 


•3535- 


17 


.3656-07 


.9451+00 


10 


.00 


X(l) = 


•2712-01 


.3474+02 


•3131- 


18 


•1256-07 


.9451+00 



END OF SIMULATION 
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APPENDIX II 



^XQT .COFLO . LCTKA - VOLTERRA EQUATIONS 



COMPARTMENT - FLOW MODEL SIMULATION PROGRAM 

• NAMES 
XlspOPi 
X2=P0P2 
•FLOWS 

F(et 1 )=Ri#pOPl*( 1 .-POP! /Kl-ALP HA* POP? /K 2) 

F{ 0t2)=R2tpOP2*( T.-P0P2/K2 -SETA •POP! /Kl ) 

'►INITIALIZE 

P0P1=13. 

H0P2=92. 

Kl =52. 

K2=120. 

ALPHA=-2. 

3ETA=.€ 

Rl=.5 

R2 = .6 

*T IMER 

ST0P=15. 

0ELTA=1. 

• RUN 

TABLE OF INCLUDED FLOWS 

0 1 2 
0 XX 
1 
2 



READY 

T= .00 X( 

T= 1.00 X( 

T= 2.00 X( 

T= 3.00 X( 

T= 4.00 X( 

T= 5-00 X( 

T= 6.00 X( 

T= 7*00 X( 

T= 8.00 X( 

T= 9.00 X( 

T= 10.00 X( 

T= n.00 X( 

T= 12.00 X( 

T= 13.00 X( 

T= 14.00 X( 

T*= 15.00 X( 



END OF PREPARATION PHASE 



.1000*02 


•9000*02 


•3263*02 


.3404*02 


.6933*02 


•6712*02 


•8958*02 


•4818*32 


.9007*02 


•3563*02 


.8494*02 


•2855+02 


.7999*02 


•2455^02 


•7626*02 


•2218*02 


•7365*02 


.2073*02 


•7180*02 


•1932*02 


•7065*02 


•1924+02 


•6982*02 


•1886+02 


•6927*02 


.1862+02 


•6890*02 


•1846+02 


•6865*02 


•1836+02 


•6849*02 


.1829+02 



END OF SIMULATION 



